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ABSTRACT 


Many proposed space reactor designs employ heat pipes as a means of convey- 
ing heat. Previous researchers have been concerned with steady state operation, 
but the transient operation is of interest in space reactor applications due to 
the necessity of remote startup and shutdown. A model is being developed to study 
the dynamic behavior of high temperature heat pipes during startup, shutdown and 
normal operation under space environments. Model development and preliminary re- 
sults for a hypothetical design of the system are presented. 


INTRODUCTION 


HEAT PIPE DESCRIPTION 


A heat pipe is an effective means of transferring heat from one location to 
another without pumps or moving parts. It consists of a closed system filled with 
a working fluid. One section of the pipe is placed in a heat source, causing the 
working fluid to evaporate and the resultant vapor to expand into the remainder 
of the heat pipe. A different section of the heat pipe is placed in a heat sink, 
causing the working fluid to condense. The liquid is then conveyed back to the 
evaporator section by capillary action through a wick structure (see figure 1). 

The primary heat transfer is due to the latent heat of the working fluid, which 
results in a nearly isothermal system. The large heat transfer with the accomp- 
anying small change in temperature makes the heat pipe equivalent to a structure 
with a very high thermal conductivity. 


DESIGN CONSIDERATIONS FOR SPACE REACTOR USE 


Many of the current space based reactor designs for both civilian and mili- 
tary applications employ heat pipes as a means of conveying heat (ref. 1). In 
these designs, thermal radiation is the principal means for rejecting waste heat 
from the reactor system, making it desirable to operate at high temperatures. 


*The research was performed under appointment to the Nuclear Engineering, Health 
Physics, and Radioactive Waste Management Fellowships program administered by Oak 
Ridge Associated Universities for the U. S. Department of Energy. 
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Lithium is generally the working fluid of choice as it undergoes a liquid-vapor 
transformation at the preferred operating temperature. There are, however, prob- 
lems inherent to the choice of a liquid metal as the working fluid. The reactor 
assembly will have to be launched into orbit, presumably in a cold state with the 
lithium solid. Therefore, the conditions under which the heat pipe will be self- 
priming from a frozen state are of great interest. A similar problem is that of 
restarting the heat pipe after it has been shut down and allowed to solidify; the 
difference between this problem and the former one being the initial distribution 
of the solid lithium. An additional concern, particularly in military applica- 
tions, is the ability of the system to handle extreme bursts of heat input con- 
centrated at a particular location such as might be expected if the satellite 
were the subject of a military attack. 


PURPOSE OF THIS WORK 


The steady state behavior of heat pipes has been studied by other researchers 
(refs. 2,3,4). However, the nature of the previously mentioned processes of re- 
mote startup, restart, and reaction to threats necessitates an accurate, detailed 
transient model of the heat pipe operation. This paper covers the development of 
a model of the vapor core region of the heat pipe which is part of a larger model 
of the entire heat pipe thermal response. Other transient heat pipe modeling has 
been done, including a thermal model of a low temperature heat pipe by Chang and 
Colwell (ref. 5) which does not include the hydraulic behavior modeled in this 
paper. Another notable effort involves a modification of the ATHENA code to model 
an entire space reactor system (ref. 6). This paper differs from the ATHENA code 
modification in its modeling of diffusion, which the ATHENA code does not treat. 

In addition, future plans for the model presented in this paper include the 
modeling of startup from a frozen state. 


MODEL DEVELOPMENT 


CORE MODEL 


The vapor core is modeled using the area averaged Navier-Stokes equations in 
one dimension, which take into account the effects of mass, energy and momentum 
transfer. The core model is single phase (gaseous), but contains two components: 
lithium gas and a noncondensible vapor. The differential form of the equation 
set used is: 


Mixture Continuity Equation 
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Noncondensible Continuity Equation 
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Mixture Internal Energy Equation 


li pu) m + ±_ 

3t 3z 


[(Pu)m v m 3 


= -P Urn 
m 3z 


+ Th; 


pc 


"(IxTrt) + f-[p,, (D N h N -D ghg l!ii!L] 

3 z 


Mixture Momentum Equation 
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State Equation 
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This equation set presupposes the thermal equilibrium of the lithium gas and 
the noncondensible, making a second internal energy equation unnecessary. The 
gases are assumed to have the same convective velocity, eliminating the need for 
a second momentum equation. However, differing diffusive velocities are allowed 
with the inclusion of concentration driven mass diffusion (ref. 7). Evaporation 
and condensation provide mass and energy transfer to and from the liquid wick 
section. Condensation is assumed to occur only at the liquid-vapor interface, 
and not by liquid droplet formation in the vapor core. The evaporation rate is 
given by the following relation: 
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Convective energy loss to the wick surface is included, but thermal conduction is 
not included in the model at present. The body forces are assumed zero for a 
space based system. 

This system of coupled nonlinear partial differential equations is solved 
using a finite difference method. First, the equations are differenced spatially, 
using a staggered mesh in which the fluid properties are evaluated at the cell 
centers and velocities are evaluated at cell boundaries. When it is necessary to 
evaluate a fluid property at a cell boundary, the variable is "donored” (symbol- 
ized by a dot over the variable) by setting it equal to the closest "upstream" 
cell center value, where "upstream" is determined by the velocity at that loca- 
tion. For example, the donored density is defined as: 
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The equations are differenced temporally using a semi -implicit discretization 
which maintains a linear system in the new time variables while adding increased 
stability over a fully explicit method. In the semi-implicit method, the velo- 
cities in the convective terms and the pressures in the momentum equation are 
evaluated at the new time step. Finally, the state equation is differenced 
separately by expanding it about the present time step in a Taylor's series, 
truncating the series after the linear terms and evaluating the truncated series 
at the new time step. The final discretized system is: 
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Noncondensible Continuity Equation 


N + l N 
P N,i P N,i 
At 


. N N+l| 
P M v 
N m 


. N N+l| 
. 1 P * 
i+^ m m 


i-12 


r N n 

1 p D 
= 1 L m 


I i+h 


Az . 

l 


Az , 


I /"N N > 

\ Az . , ^ 

i-% 



Mixture Internal Energy Equation 
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Mixture Momentum Equation 
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State Equation 
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These equations, along with a similar momentum equation evaluated at i- 7 , 
constitute a linear algebraic system of six equations in eight variables 
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Algebraic manipulation is used to reduce this to a single equation in three 
variables 
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This single equation is then applied at every node in the physical system to be 
modeled, resulting in a matrix equation with a tridiagonal structure. The 
matrix equation is easily solyed by a forward-backward sweep (Thomas algorithm) 
to yield the pressure distribution at the new time step. Back substitution 
into the discretized equation set then yields the remainder of the variables 
at the new time step. This method of discretization is similar to the methods 
used in the TRAC and RELAP computer codes (refs. 8,9). 


This finite difference method necessitates limits on the time step size in 
order to assure stability and accuracy. The Courant limit, which specifies 
that a fluid element cannot be transferred more than one cell during a time 
step, is: 
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The error inherent to truncating the state equation series expansion, referred 
to as density truncation error, is controlled by decreasing the time step when 
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the error exceeds a user specified value: 
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WALL/WICK MODEL 


In the current work, the wall and the wick of the heat pipe are modeled 
simplistically . A detailed thermohydraulic model of the wick is under develop- 
ment. The primary reasons for including a model of the wall and the wick are: 
to provide a wall temperature for calculating evaporation rates and convective 
heat losses for the core model; and to connect the core model with boundary con- 
ditions exterior to the heat pipe. Currently, the wall and the wick are modeled 
as a single entity, using a lumped parameters approach which assumes a constant 
temperature throughout the wall and the wick at any particular axial node. The 
mechanisms of heat transfer available to the wall/wick are: exterior heat input, 

exterior radiative heat loss, convection to the core gas and evaporation at the 
core-wick interface. The heat balance equation for the wall/wick is then: 
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After the core model is updated to the new time step, the wall temperature is 
up-dated by implicitly differencing this equation and solving the resultant 
nonlinear equation for the wall temperature by a Newton-Raphson rootfinding tech 
nique at every axial location. Exterior boundary conditions vary along the 
length of the heat pipe. In the evaporator section, a total net heat flux into 
the heat pipe is specified. In the adiabatic section of the heat pipe, all heat 
loss mechanisms are set to zero. Radiation is the only means of heat loss 
assumed in the condensing end of the heat pipe. 


RESULTS 


MODEL CONDITIONS 


The model has been applied to a 2.0m heat pipe that consists of an evapora- 
tor section of ze[0,.3], an adiabatic section for ze[.3,1.5], and a condenser 
section for ze[1.5,2]. These sections are used only to determine which boundary 
conditions to apply: heat input, zero heat flux, or heat output. A starting 

temperature of 700K is used. The heat pipe is operated at 15kW throughput for 
20s in order to heat it up and then the power is dropped to the desired 
operating condition of 1.184kW. A heat pipe of this design would normally 
operate at 15kW and 1500K, but the highest throughput the model will allow with 
an operating temperature of 1500K is 1.184kW if the only means of heat loss is 
radiative cooling. Initially, the heat pipe is filled with a noncondensible gas 
(air) at a pressure of 250Pa. 
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ANALYTICAL STEADY STATE RESULTS 


Most of the steady state heat pipe analysis has been based on Cotter's early 
work (ref. 3). These analyses have been primarily concerned with pointing out 
limits to the heat transfer capability of the heat pipe in the form of sonic, 
boiling, wicking, and entrainment limits. Of greater importance to this work is 
the application of Cotter's theory to predict pressure, temperature, and velo- 
city distributions inside the heat pipe's vapor core. Cotter's theory assumes 
that the temperature distribution is simply a constant value. The velocity dist- 
ribution, assuming constant heat addition in the evaporator and constant heat 
rejection in the condenser, linearly increases from zero to Q t /hfg in the evapor- 
ator, remains constant through the adiabatic section, and then linearly decreases 
to zero again in the condenser. The predicted steady state pressure distribution 
is more complicated and is presented in Figure 2, as calculated by the HTPIPE 
code (ref. 4). 


EXPERIMENTAL STARTUP OBSERVATIONS 


In his later work, Cotter describes three basic modes of heat pipe startup 
that have been observed experimentally (ref. 10). The first mode, a uniform 
startup, occurs when the vapor density is high, and there is continuum flow 
throughout the heat pipe. A second, frontal startup mode is observed when the 
vapor density is so low that free molecular flow effects are important. This 
mode has a continuum to free molecular flow transition region that moves down 
the length of the heat pipe and causes the progressing, frontal temperature 
distribution. The third startup mode is also frontal, caused now by the pre- 
sence of a noncondensible gas. This gas is swept out of the evaporator section 
and forms a plug that is compressed into the condenser end. The interface between 
the two gases is sharply defined and the temperature change across this front is 
more abrupt than that of the second startup mode. 


MODEL RESULTS 


The heat pipe was modeled for a period of greater than six minutes model 
time. The work was done in FORTRAN on a VAX-11/750 and required approximately 
36 hours of CPU time. Results from the first 20 seconds (startup results) and 
from the state of the model at 360 seconds (when many of the primary fluctuations 
have died out and the model is approaching steady state) are presented. 

Wall temperature increases evenly across the evaporator section during 
startup, with the rest of the heat pipe following behind (Figure 3). The step 
changes in temperature at the end of the evaporator are due to the simplistic 
wall/wick model where axial heat transfer has been neglected. The adiabatic and 
condenser wall sections are heated by condensation and convection from the central 
gas core. The mixture temperature generally follows the wall temperature (Figure 
4). However, the sharp temperature changes present in the wall are not present 
in the core because of the convective transport of the hot gas. The frontal 
startup shown in the model resembles a combination of the second and third startup 
modes observed by Cotter, which should be expected due to the presence of the 
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noncondensible and the low vapor densities. At 16 seconds, both the mixture and 
wall temperatures have assumed flat profiles. 

The lithium gas pressure increases with time as the wall temperature increases 
and condensation is inhibited (Figure 5). The lithium gas produced in the evapor- 
ator section travels down the heat pipe, pushing the noncondensible gas ahead of 
it (Figure' 6). By 16 seconds, the noncondensible has been packed into the end 
of the condenser with the lithium occupying the remainder of the heat pipe. The 
moving noncondensible-lithium interface is sharply defined, as observed by Cotter. 
The total pressure remains relatively flat throughout the startup transient. 

The mixture velocity exhibits a maximum during the startup of the heat pipe 
(Figure 7), with the distribution peaking at the end of the evaporator section. 

The magnitude of the peak increases quickly to its maximum value at 4 seconds, 
and then decreases slowly as lithium gas is transported further down the heat 
pipe before condensing. The maximum value of 700 m/s is much larger than steady 
state values of around 10 m/s. 

The evaporation rate during startup increases rapidly to a rather constant 
value in the evaporator section (Figure 8). Note that the 1-second line is ap- 
proximately at zero: an evaporation rate above that line indicates net evapora- 

tion, while an evaporation rate below that line indicates condensation. Due to 
the relatively cold wall temperature beyond the evaporator section, the lithium 
gas begins condensing as soon as it leaves the evaporator, with the effective 
condensing length of the heat pipe increasing with time. Condensation increases 
when the gas reaches the condenser section, which is radiatively cooled. 

After 360 seconds of model time both the wall and the mixture temperature 
show a drop of less than 20 degrees from the evaporator to the condenser (Figures 
9,10), not including the plug of noncondensible gas in the condenser end. The 
temperature variations in this plug of gas are due to relatively poor low velocity 
heat transport of the noncondensible as compared to condensation, and the lack 
of conduction in the core model which would reduce the 40 degree temperature 
gradient. The pressure distributions show that the noncondensible has been com- 
pacted into the condenser end (Figure 11). Mixture pressure is very flat, but 
does have some variation which has been magnified greatly in Figure 12. This 
distribution has the same general shape as the steady state distribution calcu- 
lated from the HTPIPE code (Figure 2). 

The mixture velocity distribution at 360 seconds consists of three linear 
sections (Figure 13). The evaporator section is linearly increasing, which 
coincides with Cotter's theory on steady state heat pipe behavior. The condenser 
section is linearly decreasing in velocity, which also agrees with Cotter's 
theory. However, the velocity profile in the adiabatic section is not flat, as 
Cotter's theory predicts, but rather decreases linearly with a different slope 
from the condenser region. The reason for this behavior can be seen in the evap- 
oration rate (Figure 14). The evaporation rate is constant in both the evaporator 
and condenser section, as Cotter assumes, but has a nonzero value in the adiabatic 
section. The nearly constant condensation rate in the adiabatic section causes 
the velocity in that section to resemble the linearly decreasing velocity of the 
condenser. Condensation in the adiabatic section indicates that the system has 
not fully achieved a steady state condition. 

Diffusion was taken out of the model in order to determine its effect on heat 
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pipe performance. Figures 15 through 20 (in which diffusion has been neglected) 
correspond to Figures 3 through 8. There is little difference between the sets 
of plots, indicating that diffusion does not affect the heat pipe performance 
during this transient. 


SUMMARY 


A transient model of the vapor core of high temperature heat pipes has been 
developed. Model trends show good agreement with both analytical steady state 
calculations and experimentally observed startup behavior. Diffusion does not 
have a large effect on heat pipe performance for this particular transient. In- 
vestigation into the transient conditions during which diffusion has a more pro- 
minent effect is planned. Future modeling efforts will include a thermohydraulic 
model of the wick region with provisions for startup from a solid state and con- 
duction in the vapor core. 

SYMBOLS 


A s - surface area of the liquid-vapor interface (m 2 ) 

Cp - specific heat (J/kg»K) 

D - mass diffusion coefficient (m 2 /s) 

- hydraulic diameter (m) 

E ma x " maximum tolerable error 

f - friction factor 

g - acceleration due to gravity (m/s 2 ) 

H - heat transfer coefficient (J/m 2 *K) 

h - enthalpy (J/kg) 

hfg - latent heat of evaporation (J/kg) 

hP c - enthalpy of the lithium gas upon changing phase (J/kg) 

M - molecular weight (kg/kg*mole) 

m — mass (kg) 

P - pressure (Pa) 

Q t - total heat throughput (W) 

- heat flux incident on the heat pipe (W/m 2 ) 

R - gas constant ( J/kg*mole *K) 

r - inner wick radius (m) 

T - temperature (K) 

t - time (s) 

u - internal energy (J/kg) 

V - velocity (m/s) 

X - mass fraction 

z - axial distance (m) 

Greek Letters: 

T - evaporation rate (kg/m 3 *s) 

A - difference between two values, e.g. AXi=X^ + l /2 - X^_]/ 2 , AXf + l/ 2 = Xi+l~Xi 
e - emissivity 

p - density (kg/m 3 ) 
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a - Stef an-Boltzmann constant (J/im 2 *s*K-^) 

Subscripts : 

g - lithium gas 

i - spatial node 

m - mixture 

N - noncondensible gas 

sat - vapor-wick saturation value 

w - wick/wall 

Superscript : 

N - time step 
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Figure 2. Steady state pressure profile (HTPIPE) 



Figure 3. Wall temperature (startup) 



Figure 4. Mixture temperature (startup) 
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Figure 8. Evaporation rate (startup) 



Figure 9. Wall temperature (360 seconds) 



Figure 10. Mixture temperature (360 seconds) 
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Figure 1 1 . Pressures (360 seconds) 
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Figure 12. Mixture pressure (360 seconds) 
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Figure 13. Mixture velocity (360 seconds) 
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Figure 16. Mixture temperature (startup, no diff) 
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Figure 17. Lithium gas press, (startup, no diff) 



Figure 18. Noncondensible press, (startup, no diff) 



Figure 19. Mixture velocity (startup, no diff) 
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Figure 20. Evaporation rate (startup, no diff) 
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